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Optimal single electron sources emit regular streams of particles, displaying no low frequency 
charge current noise. Due to the wavepacket nature of the emitted particles, the energy is however 
fluctuating, giving rise to heat current noise. We investigate theoretically this quantum source of 
heat noise for an emitter coupled to an electronic probe in the hot-electron regime. The distribution 
of temperature and potential fluctuations induced in the probe is shown to provide direct information 
on the single particle wavefunction properties and display strong non-classical features. 
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Recent years have witnessed a surge of interest in on- 
demand sources for single particles in mesoscopic and 
nanoscale systems. This interest was motivated by the 
experimental progress [H-0| on fast, accurate single parti- 
cle emitters, with operation frequencies reaching the GHz 
regime. Fast and accurate emitters are key elements in 
the efforts to obtain a quantum standard for the Ampere 
Q. In addition to metrological applications, coherent on- 
demand sources, emitting regular streams of single par- 
ticle wavepackets, are of great, fundamental importance. 
As recently demonstrated [f|, sources implemented with 
edge states in the quantum Hall regime P, [|| open up 
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for quantum coherent few-electron experiments 
well as put in prospect quantum information processing 
HI- 13] with clocked single and entangled two-particle 
sources. Large efforts have also been put into characteris- 
ing the properties of on-demand sources via the electrical 
current and its fluctuations 

Although the low-frequency charge emission of ideal 
on-demand sources is noiseless, the emitted heat fluctu- 
ates (2^ |. These fluctuations are ubiquitous for quantum 
coherent sources; particles emitted during a time shorter 
than the drive period T have an uncertainty in energy 
larger than H/T. Acting as emitters of quantum heat 
fluctuations, coherent on-demand sources comprise ideal 
components for tests of heat fluctuation relations [28 - 31 1 
in the quantum regime I32B or for investigating the statis- 
tics of temperature [33l |34| or heat transfer [3^| fluctua- 
tions in mesoscopic systems. In addition, the large versa- 
tility of system parameters and pulse protocols p, Q for 
on-demand sources allows for a tailoring of the spectral 
properties of the emitted wavepackets. 

In this work we provide a compelling illustration of the 
fluctuation properties of system consisting of a generic co- 
herent on-demand source coupled to a hot-electron probe, 
see Fig. [TJ It is shown that the temperature and poten- 
tial fluctuations induced at the probe, besides fundamen- 
tal constants, depend only on the source frequency and 
the spectral properties of the emitted wavepackets. For a 
wide range of parameters, the quantum fluctuations are 



found to dominate over the classical ones. In addition, 
the full distribution of the fluctuations reveals a direct 
proportionality between the cumulants of the marginal 
temperature and potential distributions, allowing for an 
experimental investigation of the temperature fluctuation 
via correlators of the potential fluctuations. 
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FIG. f : a) Schematic of an on-demand source injecting single 
particle wavepackets into an electronic probe, via the lower 
edge state of a conductor in the quantum Hall regime. The 
probe is in the hot-electron regime, with a floating electron 
temperature T p (t) and chemical potential fi p (t). Particles 
emitted from the probe flow along the upper edge into an elec- 
tronic reservoir electrically grounded and kept at zero temper- 
ature, b) Noiseless train of wavepackets, emitted from the sin- 
gle particle source with a frequency u — 2tt/T. c) Probability 
distribution of energy p(e) of the wavepacket, with average (e) 
and width Ae shown. 

We first discuss the energy emission properties of an 
isolated, optimal on-demand source, emitting a train of 
single particle wavepackets \ty), equally spaced T = 2-k/uj 
in time, see Fig. [1] The wavepackets, emitted on top of 
a filled Fermi sea, are superpositions of states at different 
energies, 



I*) = 



dec(e)tf(e)\0), 
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where w(e) creates a particle at energy e > 0, |0) denotes 
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the filled Fermi sea and c(e) an amplitude normalised as 
Jdep(e) = l, withp(e)=Jc(e)| 2 ._ 

Although energy emission for individual electrons can 
be accessed via charge counting in weakly tunnel coupled 

systems m, mi, [Hi , for on-demand sources operating in 

the GHz regime single-shot energy detection is presently 
not possible. Instead one has to consider schemes where 
the effect of collecting a large number N ^> 1 of emit- 
ted electrons, with a fluctuating total energy E, become 
measurable. For an on-demand source characterized by 
p(e), the statistical distribution P(E) of the total energy 
can conveniently be written 



P(E) = J d 



i£E+NF{Z) 



F(0 = In 



(2) 



where NF(£) is the cumulant generating function for 
P(E) and N = to/T is the number of particles with 
to 3> T the measurement time. The different cumulants 
of P(E) are obtained by successive derivatives of F(£) 
with respect to £, giving for the average and the width 

(E)=N(e), (AE) 2 =N((e 2 )-(e) 2 )^N(Ae) 2 (3) 

where (..) = J de..p(e). Importantly, the direct relation 
between the statistics of E and p(e) of the individual 
wavepackets depends crucially on the ideal operation of 
the source. Irregular wavepacket emission or scattering 
in space or energy between emission and detection will 
make P(E) dependent on other factors. We stress that 
on-demand sources which emit subsequent, identical elec- 
tron and hole wavepackets [H, (and hence no net charge) 
have the same heat emission properties as pure electron 
sources with twice the drive frequency oj. 

To access the heat fluctuation properties of the source 
we consider an on-demand source coupled to a probe 
in the hot-electron regime, an electrically and thermally 
floating terminal, see Fig. [TJ The source-probe setup is 
implemented in a conductor in the quantum Hall regime, 
allowing to minimise scattering, elastically or inelasti- 
cally, between particle emission and collection in the 
probe. Particles emitted from the source propagate to 
the probe along the lower edge state. From the probe, 
emitted particles follow the upper edge and are collected 
in a grounded electronic reservoir kept at low (here zero) 
temperature. 

In the hot-electron probe, injected particles thermal- 
ize rapidly, via electron-electron interactions, on the time 
scale T e _ e . This time is much shorter than the typical 
time Td a particle spends inside the probe before being 
reemitted. However, the energy exchange with the lat- 
tice phonons takes place on the time scale T e - P h, much 
longer than Td- The electron distribution in the probe 
is then in a quasi-equilibrium state, characterized by a 
chemical potential and temperature. In order to prevent 
charge and energy pile-up in the floating probe, both the 
chemical potential fi p (t) and temperature T p {t) develop 



fluctuations in time. The potential fluctuations can be 
detected by present day electrical measurements. Impor- 
tantly, as we now show, the potential fluctuations also 
provide direct information on the quantum heat fluctua- 
tion of the source, via p(e) of the individual particles. 

To present a clear and compelling picture we anal- 
yse the temperature and potential fluctuations within 
a Boltzmann-Langevin approach. We write T p (t) = 
T p + ST p (t) and n P (t) = p p + 5/j, p (t) with T p ,pL p av- 
erage quantities and ST p (t),5fi p (t) fluctuating Langevin 
terms. The statistics of 6T p (t) and Sp p (t) is determined 
from underlying quantum properties, as discussed be- 
low. First we consider the average quantities T p and 
pp. The starting point is the operator for charge current 
at the probe [Hi, J|(i) = / dede' e i{ - e - e 'W h i p (e, e') where 

i p (e,e') = [e/h)p{e)b{e') - a){e)a{e')] with fr(e)[tf(e)] 
creating particles incident on [emitted from] the probe. 
By analogy we write the operator for the heat current 
f£(t) = J dede'e^'^^Ke + e')/2 - fj, ]i p (e,e'), where 
fiQ = is the chemical potential of the reservoir. Taking 
the quantum average with respect to the emitted source 
state, Eq. {1} , and the state of the probe, and averaging 
over a time much longer than the period T, we arrive at 
the dc component of charge and heat currents 
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(4) 



where go = e 2 /h is the (single spin) conductance quan- 
tum and Iq = (irks /e) /3 the Lorentz number. To ac- 
count for both types of sources discussed, we introduced 
(7 = for sources emitting no net charge and a = 1 for 
sources emitting one electron per cycle. We note that 
the first and the second terms of (I p ) and in Eq. 
(HJ are the currents emitted by the source and the probe 
respectively. The conditions for zero average charge and 
heat currents at the probe, (I p ) = and (I p ) = 0, give 
from Eq. (gj) 



/i p = aftw, T p = 



1 



goloT 



[2(e) - ahuj]. 



Importantly, p p and T p depend only on the source prop- 
erties uj and (e) and fundamental constants [j|. We note 
that (e) > fku/2 follows along the lines of Ref. [37| . 

Turning to the temperature and chemical poten- 
tial fluctuations [38[, the quantities of primary ex- 
perimental interest are the low frequency correlators 
((Sfip) 2 ) = (1/to) / dtdt' (Sp p (t)5fi p (t')) and equivalently 
for ((ST p ) 2 ), with the measurement time to ^> T. We 
first point out that the total fluctuations of charge and 
heat currents Alp and Alp are made up by bare fluctu- 
ations Sip and dip and fluctuations due to the varying 
temperature and voltage of the probe, dr p {I p )8T p and 
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suppressing for shortness the time dependence of the fluc- 
tuations. From Eqs. (j4} and ([6]), taking into account the 
conservation of charge and heat current fluctuations, i.e. 
Alp = and AI p — 0, we can express Sfj, p and 8T p in 
terms of the bare charge and heat fluctuations as 
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8ii p = ~5i;, ST p = —^ r (SI^-^5I p A. (7) 

The correlators ((8^, p ) 2 ) and ((8T p ) 2 ) can be thus 
expressed in terms of low frequency correlators 
of bare charge and heat fluctuations (51*81%) = 
(1/to) / dtdt'(8I^(t)Siy{t')). The correlator of the 
Langevin terms (8I p (t)8I%(t')) is evaluated by taking the 
quantum average of the corresponding correlator of cur- 
rent operators I p , I p following Ref. [36(. We arrive at 



((8fi p ) 2 ) = hk b T p , 



(8) 
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The potential fluctuations ((8fi p ) 2 ) are proportional to 
the average temperature, typical for equilibrium systems 
[38l |. In contrast, the temperature fluctuations ((ST p ) 2 ) 
are a sum of two physically distinct terms. The first term, 
the classical fluctuations, is proportional to T p and results 
from the finite temperature of the probe and would be 
present even if the injected particles had a well defined 
energy, i.e. Ae = 0. The second term, quantum fluc- 
tuations, is proportional to (Ae) 2 and is a direct result 
of the uncertainty of the energy of the injected particle. 
Importantly, for a broad range of drive frequencies u) and 
wavepacket mean energies (e) and widths Ae, the quan- 
tum term dominates over the classical one. We also point 
out that there are no correlation between the voltage and 
the temperature fluctuations, i.e. (8/i p 8T p ) = 0. 

In order to investigate the presence of quantum heat 
fluctuations in higher order potential correlations, and 
also to provide a complete picture of the temperature 
and potential fluctuations, we turn to the full proba- 
bility distribution. To relate to the average and fluc- 
tuation correlators above, we introduce a dimension- 
less potential /i = (1/h) J ° dtfj, p (t), and temperature 

T = (1/h) L° dtkbTp(t), fluctuating quantities integrated 
over the measurement time to . The joint probability dis- 
tribution Vt (/■*, T) can be conveniently written in terms 
of a cumulant generating function G(x, 9) as 



1 



(27T) 



d X / d0 e -* r -*M+G(x,9) ) (9) 



with x ancl $ counting fields for [i and T respectively. 
From G(x, 9) the low frequency cumulants are then, by 



construction, obtained from successive derivatives with 
respect to the counting fields, giving to(8T p 8 /j,™ 1 ) — 

(-ihr+ m k^ n d2d™G(x,6)\ x ,g =0 . 

To determine G(x, 9) we first spell out the relations 
between the time scales in the problem. The potential 
jUp(i) fluctuates on the time scale given by the RC-time, 
trci while the temperature T p (t) typically fluctuates on 
the time scale of the dwell-time in the probe, r^. We 
assume that the system is in the limit r e _ e -C Trc, id *C 
Te-pft. Moreover, we consider periods of the source, T, 
and measurements time such that to 3> t^, trc ^ T. On 
time scales r such that T <C t < r,;, trc the statistics of 
net transferred energy E p and charge Q p in the probe can 
be described by the source generating function rh s (X, £) 
with 



h. = 



2- 



[-iea\ + F(f)] 



(10) 



with F(£) given in Eq. $2$ and the probe generating 
function r/i p (A,^, E p ,Q p ) with I3T 



41] 



ln[l + / p (e)(e 
ln[l + /(e)(e 



-ieX— ie£ 
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-1)] 



(ii) 



where / p (e) = f p (e, fj, p ,T p ) and /(e) are the probe and 
the reservoir distribution functions (see Fig. Q} and £ 
and A are the counting fields for E p and Q p respectively. 
The energy E p and charge Q p are related to T p and fj, p 
as E p = v[/j. p /2 + (irkbT p ) 2 /6] and Q p = ven P , where 
v is the density of states in the probe. Working within 
the framework of the stochastic path integral formalism 
41, 4^, we can then express G(x,9) as a path integral 
over all configurations of E p and charge Q p during the 
measurement. In the long time limit we have 



J dQ p dE p d\d£e s( ' Q '" E '" x >t ) (12) 



33j,l4fJ,l41|, 



where S(Q P ,E P ,\,£) = t [i9k b T p /h -I 
h p (Q p , E p , A, £) + h s (X, ^)]. Similar to Refs. 
the integral in Eq. (fT2|) is solved in the saddle point ap- 
proximation. Inserting the solutions for Q pi E p , A,£ into 
, A, £) we arrive at 



S(Q P , E p , . 



G( X ,0) = N 



d[zF(z)\ 
dz 



+ a(z + ix) 



(13) 



recalling that N = t Q /T, where F(z) = F(^)\^ =z/hw 
and z is found from the relation z 2 [dF/dz + a/2] = 
-(tt 2 /6)(1 - yi— 2^) 2 with g( X , 9) = (3/7r 2 )[( iX ) 2 /2 + 
i9}. From Eq. (fl3|) we note several important things. 



First, by expanding G(x, 9) in terms of x an d 6 we see 
that the first two cumulants reproduce the results in Eqs. 
(0) and ([8]). Second, all the even chemical potential cu- 
mulants, from the first two terms in Eq. (fT3"|) . can be 
expressed in terms of the temperature cumulants as 



(2n-l)\l(k b hT((8T p ) n ) 



(14) 
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a consequence [43| of the counting fields entering via 

g(x,d)- 

The full distribution Vt {fx, T) can be found (to expo- 
nential accuracy) by solving the integral in Eq. © in 
the saddle point approximation. We obtain the compact 
expression 

InVtMT) = -iT6* + G(0, 6*) — (fJ, — H) 2 /{2T) (15) 

where p, = toflp/h = aN and the saddle point solu- 
tion 8* is found from the relation dF/dz\ z=z * + a/2 = 
-(TT 2 /6)( Tq*/N) 2 , w ith z* = N(q* - l)/{q*T) and 
q* = \Jl — i68*/n 2 . Importantly, Eq. (TT5"]) shows 
that the potential fj, displays Gaussian fluctuations, of 
width VT, around the average p, for any given tem- 
perature T. Hence, the marginal potential distribution 
Vt (^) = f dTVt (fi>TT) is symmetric around fx, albeit 
not Gaussian. Moreover, the marginal distribution for 
the temperature Vt {T) = J d[j,Vt {n,T) is given by 
lnV to (T) = -iT6* + G(0,6*). 

To illustrate our results we evaluate In "Pt (T) for 
two distinct cases. First, as a reference, we con- 
sider the generic Gaussian spectral distribution p(e) = 
l/(\/2^Ae)e-( e -< e » 2 /( 2 ( Ae ) 2 >. Taking the classical limit, 
with Ae <C (e) , we get the simple result 



lnV t0 (T) = -{-k 2 /Q)T{1-T/T)' 



(16) 



with T = takbTp/h the average value of T. The log proba- 
bility is plotted in Fig. O For small fluctuations T — T <C 
T the distribution is Gaussian while for T <C T the prob- 
ability is suppressed Vt (T) oc e~^' T ~ '/( 6T \ guaranteeing 
V to (T) -> for T -> 0. The probability for large fluc- 
tuations T > f is suppressed as T> to {T) oc e~ T?r2/6 . For 
finite but small width Ae <C (e) , the log probability in Eq. 
(TTo) is multiplied by the term l + (Ae/(e)) 2 T 2 7r 2 /(12TiV) 
for T ~ T, a small quantum correction. In contrast, for 
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FIG. 2: Normalised logarithm of the probability distribution 
Vt (T) as a function of T/T, with T the average value of T. 
The curves correspond to the narrow Gaussian wavepacket 
energy distribution (blue solid) in Eq. (|16[1 and to an expo- 
nential energy distribution with a — 1 (black dotted), 3 (red 
dash-dotted), 5 (green dashed). See text for details. 



an exponential distribution p(e) = (l/(e))e -e /^ derived 
in Ref. (l9| for adiabatic particle emission and investi- 
gated in Ref. [l^, we find the probability 



ln7> t0 (T) = - 



TT 2 T 

~6~ 



T 2 
I a 



Tq* 



T 



(17) 



where q* = (a/2 + yJaT/T + 1 + a 2 /A) /{T/T + a) and 
a = n\J (e) / (6hiu). The probability distribution is plot- 
ted in Fig. [2] for different values of a. We see that for in- 
creasing a, corresponding to slower drive ui and/or larger 
average wavepacket energies (e), the distribution gets in- 
creasingly broad and deviates strongly from the classical 
one in Eq. (fT6)) . 

In conclusion we have investigated the quantum fluctu- 
ations of the heat current emitted from a single particle 
source. We show that these quantum heat fluctuations 
can be detected via electrical potential fluctuations of a 
probe coupled to the source. For typical parameters [l| 
2iru) r~j 1GHz and (e) ~ O.lmeV we get fl p ~ 5fj,eV and 
T p = 0.2K, demonstrating the experimental feasibility of 
our proposal. 
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Swedish VR. 
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